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Abstract- The paper describes a new method for determining 
critical load for a bus with respect to its voltage collapse limit of 
an interconnected multi-bus power system. For this purpose, 
network partitioning technique is used to transform load flow 
Jocobian matrix into a two by two matrix with respect to a 
selected/target bus and a procedure has been developed to 
determine load margin of the target bus using the elements of the 
transformed two by two elements Jacobian matrix and the bus 
voltage of the target bus. The validity of the proposed method 
has been investigated for the IEEE 30 and IEEE 118 bus systems. 

Keywords- Voltage Stability; Load Margin; Voltage Collapse; 
Power System 

List of Symbols: 

N= Total number of buses in the system 

NG = Total number of generation buses in the system 

Pi = Injected active power at ith bus. 

Q i = Injected reactive power at ith bus . 

P Di = Active power demand at ith bus. 

Qdi= Reactive power demand at ith bus. 

Poi = Active power generated by ith generator. 

Qci = Reactive power generated by ith generator. 

Vi = Magnitude of voltage at ith bus . 

5i = Angle of the bus voltage at ith bus. 

cos = Load power factor of the ith load bus. 

Qj + j By = Element of Y-BUS matrix at ith row and jth 

column. 

I. INTRODUCTION 

The voltage instability of an interconnected power system 
is a major concern for the power system operator and planner. 
It has been observed that voltage magnitudes in general, do 
not give a good indication of proximity to voltage stability 
limit In recent literature several voltage stability and 
voltage collapse prediction methods have been presented. 
Some of the important ones are: voltage collapse index based 
on a normal load flow solution [2 ' 3,4] , voltage collapse index 
based on closely located power-flow solution pairs [5] , 
voltage collapse index based on sensitivity analysis [6] , and 
minimum singular value of Newton - Raphson power flow 
Jacobian matrix [7, 8] . Continuation power flow analysis is 
based on locally parametrized continuation technique; it aims 



at avoiding the singularity of the Jacobian by slightly 
reformulating the power flow equations [9] . Network 
partitioning technique has been employed for investigating 
voltage stability condition of a load bus [10] . A simple, 
computationally very fast local voltage-stability index has 
been proposed using Tellegen's theorem It is easy to 
implement in the wide-area monitoring and control center or 
locally in a numerical relay. A new node voltage stability 
index called the equivalent node voltage collapse index 
(ENVCI), which is based on ESM and uses only local voltage 
phasors, is presented . Information/knowledge of load 
margin of a bus with respect to its voltage collapse limit 
constitutes an important criterion for load pick-up step during 
power system restoration planning and operation. The bus 
having highest load margin appears to be the best choice [13] 
for load pick-up and magnitude of load to be picked-up must 
be less than the load margin for a bus. 

A method has been proposed in reference [14] to 
determine load margin of a target bus by transforming the 
load flow jacobian matrix of a multi-bus system into a two by 
two matrix with respect to the target bus. But, the 
transformation step uses sensitivity relation which depends 
on decoupled basis (simplification). In addition to this, the 
perdition step uses an empirical formula to normalize the 
effect of linear representation of the sensitivity relations. As a 
result, the perdition defers from actual value when perdition 
is made at a point which is not near to the proximity of 
voltage collapse. This is the major disadvantage of the work 
and therefore, it is not reliable for practical use. This paper 
proposes an algorithm for determining the load margin of a 
selected/target bus with respect to its voltage collapse limit 
using the two by two matrix created by network partitioning 
technique. The use network partitioning technique eliminates 
the simplification adopted in the reference [14] while 
transforming the load flow jacobian matrix of a multi-bus 
system into a two by two matrix with respect to the target bus. 
Also, the paper introduces a full algorithm for the 
determination of load margin irrespective of working point of 
the power system Two factors namely RF (Reduction Factor) 
and DF(Distance Factor) have been introduced for this 
propose. This allows the incorporation of a correction load 
flow step after perdition step. 
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E. VOLTAGE STABILITY CRITERIA : NETWORK 
PARTITIONING TECHNIQUE 

The determinant of load flow Jocobian matrix reduces as 
an interconnected power system approaches the proximity of 
voltage collapse and it becomes zero when the system arrives 
at the point of voltage collapse. Network partitioning 
techniques [10] is used to transform a load flow Jocobian 
matrix into a two by two matrix with respect to a 
selected/target bus. As such, when this two by two matrix 
becomes zero, the system voltage collapse occurs. 

The liberalized system of load flow equation for an 
interconnected power system is represented as : 



(i) 

, where 

E, F, G and H corresponds to the load flow Jocobian of the 
system Rearranging items by putting the equations dealt with 
the bus under analysis (selected/target kth bus) to the bottom of 
the Jacobian matrix, we have, 
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, where - 

sub matrices A, B, C, D is originated from the system Jocobin 
matrix. It is assumed that the active power and reactive power 
of all the load (or generator) keep constant except for 
selected/target kth bus, i.e., AP ? =0, AQ f =0, the matrix 
represented by Equation (2) is simplified as [10] : 
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, where 
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m. DETERMINATION OF LOAD MARGIN 

Equation (3) can be represented as two bus system with the 
target kth load bus, Y-Bus elements and an equivalent source 
as depicted in Figure 1. 
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Figure 1: Equivalent two bus system with target kth bus 
Y-Bus elements and an equivalent source. 

The expression for real and reactive power injections at the 
target kth bus of the equivalent two bus system can be 
expressed as: 

, where- Gkk, Bkk, GkG and BkGare the elements of admittance 
matrix [Y] for the equivalent two bus system 
The elements of Jacobian matrix for the equivalent two bus 
system can be expressed as 



m 
m 

dQk 

m 



VcVkiGacosW + BKMh)) (10) 



At the point of voltage collapse the determinant of the 
Jacobian matrix becomes zero, i.e., 

The change in determinant value of the two bus equivalent 
system of the multibus system with respect to change in Vk and 
5k can be expressed as 

oik <'l j. 



f ftfo 

[ 3ST wt 
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(13) 
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Again, the expression for J D in terms of ^tffe 5 Wfc * 

and is as follows 
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Therefore, can be expressed as 



Pm + jQm 
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06 k 
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(15) 
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The terms 3<5jt 1 Mfc * 85* and are derived 

utilizing Equations (8) to (11). It is to be ensured that these 
terms contain only the element of the modified two by two 
element Jacobian matrix and the bus voltage of the target bus, 
because,- they are the known values for the transformed two by 
two systems. 
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Now, applying these terms of Equations (16) 
Equation (15), we have 
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can be expressed as 

Mir 



Again, utilizing Equations (8) to (1 1), we have 
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Processing Equation (9) and Equation (10), Gkk can be 
expressed as 



Oil 1 <Xj*> 



(26) 



Processing Equation (8) and Equation (11), B kk can be 
expressed as 
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Now, applying values of Equations (22), (23), (24), (25), 
(26) and (27) in Equation (21) we have 
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The terms 35^ ' 31 * " <M fe and £>V(, of Equations 
(20) and (28) are the elements of the modified Jacobian 
matrix of a multibus system represented by Equation (3) and 
Vk is the bus voltage of the kth target bus. Now, to relate 
change in real power injection at target kth bus to the change 
in determinant of the modified Jacobian matrix of a multibus 
system, variables A8k and AVk of Equation (13) are replaced 
by AP k and AQ k as follows : 



1 I 



a At 



Mi JET 

= /lAft + /a/fln^A/\ <fbr load pf same 
K/i + /i*™<MA/WAJ\ {29) 

The system collapse occurs when J D becomes zero. Now, 
if J D ° is the determinant value corresponding to the defined 
operating condition of the system, then the required change in 
AJ D of Equation(29) is AJ D = - Jd° and corresponding 
change in injection at kth bus will be 



— 1° 



(30) 



Therefore, load margin at kth load bus, i.e., the additional 
load that can be supplied to the kth load bus to push it to the 
proximity of voltage collapse is:- 



pjnw = _ APL = ^ 



f 



(31) 



Therefore, predicted critical load for kth load bus at the 
point of voltage collapse is:- 

But, P Dk margin is determined using linear relation between 
AJ D and AP k . As such, load margin for the kth load bus 
p^maigm would be more than the actual load margin of the 
bus. Further, for wide change in AJ D the predicted load value 
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will be considerably high compared to that of actual critical 
load value at the point of voltage collapse. Therefore, it is 
required to confirm the actual load margin or critical load of a 
load bus through an iterative load flow analysis using the 
load margin determined by the Equation (31). 

To ensure convergence of the load flow analysis in the 
proposed iterative procedure, the predicted critical load value 
p^margin | s norma ii ze( i ^ sucn a wa y that the modified 

predicted critical load value remains below the actual critical 
load of the bus. For this purpose, a reduction factor (RF) and a 
distance factor (DF) are used to normalize the predicted critical 
load value for the kth load bus as follows: 



~DF ™ 



(33) 



The distance factor (DF) will be higher, when prediction is 
done for wide change of AJd, i.e., far from the proximity of 
voltage collapse limit. Thus, it will normalize the effect of over 
prediction due to wide change of AJ D . Therefore, as the 
iterative procedure approaches, the proximity of voltage 
collapse limit P Dk margin will become very s mall; thus, the term 

will also approach 1. When the value of 
Jd becomes very small (say 0.0001 or less), the iterative 
procedure has to be terminated. 

IV. PROCEDURE FOR DISTRIBUTION OF NORMALIZED 
PREDICTED LOAD MARGIN TO THE GENERATORS 

In a practical system, the increase in load at a bus is 
contributed by more than one generator. Therefore, it is 
proposed to redistribute the normalized predicted load margin 

RF rymargin 
UF* Die 

of the kth target bus among the generators 
using their moment of inertia as the basis. Therefore, 
expression for the modified generation of ith generator is as 
follows : 




(or i = L .sWG (35) 



, subjected to the limiting constraints 
Pd< Pcf* for i = X...NG 



(36) 



, where- H, and P Gi max are the moment of inertia and the 
maximum limit on generation for ith generating station, 
respectively. K is the iteration count. 

V. ALGORITHM OF THE ITERATIVE PROCEDURE 

Reduction factor (RF) is used to normalize the over 
prediction of load margin (which is determined using linear 
relation between AP Dk margin and AJ D governed by Equation 



(29) with the objective of ensuring convergence of load flow 
analysis used in the proposed iterative procedure. As such, 
load at a bus must be always less than actual critical load of 
pekt 

the bus, i.e. to say that P* determined using Equation (33) 
must be less than actual critical load of the bus. In case, 
perl 

Dk determined using Equation (33) becomes slightly more 
than actual critical load of the bus (due to improper selection 
of Reduction Factor (RF), the load flow analysis of the 
iterative procedure will not converge. To take into account of 
such a situation, the proposed algorithm is equipped with a 
step after the load flow analysis. This step reloads the loads 
and generations of the system of the previous iteration values 
and reduces the reduction factor (RF) as RF = 1/1.5 and load 
flow analysis is carried out again before proceeding to the 
other portion of the algorithm This step helps in changing the 
value of reduction factor (RF) to ensure proper normalization 
(by Reduction Factor (RF)) of critical load governed by 
Equation (20). The algorithm of the proposed iterative 
procedure is as follows :- 

1. Initialize the load flow data for the system and set 
iteration count K = 0, RF = (0.3 to 0.5). 

2. Initialize the bus voltages [V] and angles [6] and conduct 
load flow analysis of the system 

3. Check for load flow convergence criteria. If load flow has 
converged, then go to Step-4. Otherwise, reduce reduction 



factor RF = 1/1.5 and reload 



pK+1 _ pK 



and 



Qff-M = qK 

**Dk **D*. Also the generating stations output with 
their Kth iteration values and go to Step-2 
pert 

4. Determine . Dk using Equation (33) and assign 

^ +1 = ^q2» +1 ='™a*& +1 

HF prnargin 

5. Distribute the additional load DF- Dk among the 
generating stations based on defined criteria subjected to 
generation limits of the generating stations as described in 
Section IV. In case of reactive power limit violation, a PV 
bus has to be changed into a PQ bus by assigning reactive 
power at its limit. 

6. Check for I -?d\ < *(= 0.0001), g0 to Step 7. If not, 
set K= K + 1 and go to Step 2. 

7. Stop. 

VI. SIMULATION, RESULTS AND DISCUSSIONS 

To verify the validity and applicability of the proposed 
method, simulations were carried out on IEEE 30 and IEEE 
118 bus systems. The proposed algorithm is used to 
determine critical load of a bus with different reduction 
factors (RF). It has been found that RF value between 0. 1 to 
0.5 ensure convergence of load flow analysis in the iterative 
process for IEEE 30 and IEEE 1 18 bus systems for any target 
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bus of the system. But, with RF = 0.3 to 0.5, the iterative 
procedure terminates with less number of iterations. During 
iteration, the additional load assigned on the target (k) bus 
was distributed among all the generating stations based on 
their load contribution (subjected to their limits) as described 
in Section 3. Continuation power flow analysis technique is 
also used to determine the critical load at the target buses 
with step size o = 0.001, because the correction load flow 
does not converge for step size greater than 0.001. 

Table 1 represents the simulation results for some of the 
buses of IEEE 30 bus system with RF = 0.5 with power factor 
0.8. Table 2 represents the simulation results for some of the 
buses of IEEE 118 bus system with RF = 0.5 with power factor 
0.8. 
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PM repesents proposed method and CLF represents 
continuation load flow. 

Tables represent the target bus selected for determination of 
critical load with respect to its voltage collapse limit, initial 
value of two by two Jacobian matrix for the target bus, initial 
value of predicted continuation parameter (AX), initial load at 
the target bus, final value of two by two Jacobian matrix, final 
value of predicted continuation parameter (at the point of 
collapse), critical load for the target bus and remarks about the 
methods . 

It is observed that J D value of two bus equivalent system of 
IEEE 30 and IEEE 118 bus systems have different initial 
values for different target buses, as such, J D value for a bus 
reflects the voltage stability characteristic of the bus. 



It has been observed that the distance factor DF = 

\ / becomes very close to 1, when the 

iterative process terminates. But, at the beginning of the 
iterative process it appears to be high depending upon 
change of AJD used for the prediction of load margin for the 
load bus. It normalizes the effect of over prediction of load 
margin due to wide change of AJD and ensures convergence 
of the load flow analysis of the system during the iterative 
process. 

The programs were executed on a P C with Pentium- 4 
processor having processor speed of 1.5 GHz and LINUX 
operating system The simulation results show that the 
proposed method requires considerably less CPU time 
compared to that of continuation power flow analysis. 
Continuation power flow analysis requires on an average 
2(Approx) times more CPU time than that of proposed 
method. 

VII. CONCLUSION 

The paper proposed an algorithm for voltage stability 
analysis of a target/selected load bus using the singularity 
condition of the load flow Jcaobian matrix. Network 
partitioning technique is used to transform load flow 
Jocobian matrix into a two by two matrix with respect to a 
selected/target bus. It is observed that determinant value of 
two bus equivalent system of IEEE 30 and IEEE 118 bus 
system have different initial values for different target buses, 
as such, it shows that the transformed two by two elements 
Jacobian matrix reflects the property/quality of the target bus. 

The algorithm proposed for the determination of critical 
load of a bus with respect to its voltage collapse limit of a 
power system works for all buses of IEEE 30 and IEEE 118 
bus systems, as such, it could be used for any interconnected 
power system The use of reduction factor (RF) and distance 
factor (DF) ensures convergence of the load flow analysis of 
the system during the proposed iterative procedure. These 
two factors effectively normalize the prediction of load 
margin, which is carried out using the linear relation between 
AJ D and AP k governed by equation (33). The simulation 
results show that the proposed method requires considerably 
less CPU time compared to that of continuation power flow 
analysis. 
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